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ABSTRACT 

NGC5471B has been suggested to contain a hypernova remnant because of 
its extraordinarily bright X-ray emission. To assess its true nature, we have 
obtained high-resolution images in continuum bands and nebular lines with the 
Hubble Space Telescope, and high-dispersion long-slit spectra with the Kitt Peak 
National Observatory 4-m echelle spectrograph. The images reveal three super- 
nova remnant (SNR) candidates in the giant Hll region NGC5471, with the 
brightest one being the 77 x 60 pc shell in NGC5471B. The Ha velocity pro- 
file of NGC 5471B can be decomposed into a narrow component (FWHM = 41 
km s _1 ) from the background Hll region and a broad component (FWHM = 148 
km s _1 ) from the SNR shell. Using the brightness ratio of the broad to narrow 
components and the Ha flux measured from the WFPC2 Ha image, we derive 
an Ha luminosity of (1.4 ± 0.1) x 10 39 ergs s" 1 for the SNR shell. The [S n] 
A6716/A6731 doublet ratio of the broad velocity component is used to derive an 
electron density of ~ 700 cm -3 in the SNR shell. The mass of the SNR shell is 
thus 4600 ± 500 M . With a ~ 330 km s _1 expansion velocity implied by the 
extreme velocity extent of the broad component, the kinetic energy of the SNR 
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shell is determined to be 5 x 10 51 ergs. This requires an explosion energy greater 
than 10 52 ergs, which can be provided by one hypernova or multiple supernovae. 
Comparing to SNRs in nearby active star formation regions, the SNR shell in 
NGC5471B appears truly unique and energetic. We conclude that the optical 
observations support the existence of a hypernova remnant in NGC5471B. 

Subject headings: galaxies: individual (M101) - Hll regions - ISM: bubbles - 
ISM: kinematics and dynamics - supernova remnants - X-rays: ISM 

1. Introduction 

Two supernova remnants (SNRs), MF83 (Matonick & Fesen 1997) and NGC5471B 
(Skillman 1985), in the giant spiral galaxy M101 have been recently suggested to be can- 
didates for "hypernova remnants" (Wang 1999), because their explosion energies appear to 
be 1-2 orders of magnitude higher than the canonical 10 51 ergs explosion energy for normal 
SNRs (Jones et al. 1998). Such high explosion energies are rivaled only by the "hypernovae" 
proposed by Paczyhski (1998) and Fryer & Woosley (1998) to explain the gamma-ray bursts 
(GRBs). As the explosion energies of MF83 and NGC5471B were derived from their X-ray 
luminosities measured from ROSAT observations which did not resolve these sources spa- 
tially, it is of interest to investigate the nature of these two hypernova remnant candidates 
at optical wavelengths with high angular and spectral resolutions. 

The analysis of optical and X-ray observations of MF 83 has been reported by Lai et al. 
(2001, hereafter Paper I). It is shown that MF83 is a superbubble around an OB association 
in which multiple supernova explosions are likely to have taken place. Optical observations 
do not provide compelling evidence for a single, super-energetic supernova explosion. Fur- 
thermore, Chandra X-ray observation of M101 shows that MF83 is dominated by a variable 
point source (Snowden et al. 2001). Therefore, it is unlikely that MF83 is a hypernova 
remnant. 

The hypernova remnant candidate NGC 5471B was initially identified and confirmed as 
an SNR in the B-component of the giant H II region NGC 5471 (hence the name NGC 5471B) 
by Skillman (1985) based on its nonthermal radio emission and high [S ll]/Ha ratio. The 
shocked, high-velocity gas in NGC5471B was detected by Chu & Kennicutt (1986), who 
concluded that the SNR was supermassive. The X-ray emission from the giant H II region 
NGC 5471 was reported by Williams & Chu (1995), but the association of X-ray emission 
with the SNR per se was made by Wang (1999). 

NGC 5471B is an unusual SNR in an active star formation region. To study the giant H II 
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region NGC5471 and the SNR NGC5471B, we have obtained Hubble Space Telescope (HST) 
Wide Field and Planetary Camera 2 (WFPC2) images of NGC 5471 in emission lines and 
continuum bands. These images are used to search for SNR candidates in NGC 5471, and to 
examine the physical structure and stellar/interstellar environment of the SNR NGC5471B. 
We have also obtained new echelle spectra of NGC5471B with the 4 m telescope at Kitt 
Peak National Observatory (KPNO) in order to kinematically separate the SNR from the 
background Hn region. This paper reports these observations (§2), the SNR candidates 
identified in NGC 5471 (§3), a detailed analysis of the physical properties of NGC5471B 
(§4), and a discussion of the nature of NGC5471B as a hypernova remnant (§5). 

2. Observations & Data Reduction 

2.1. HST WFPC2 Images 

The HST WFPC2 images of NGC 5471 were obtained on 1997 November 1 for the Cycle 
6 program GO-6928. The observations were made through two emission-line filters, F656N 
(Ha) and F673N ([S n]), and two continuum filters, F547M (Stromgren y) and F675W 
(WFPC2 R) (Biretta et al. 1996). The journal of observations is given in Table I. 

The calibrated WFPC2 images, produced by the standard HST pipeline processes, were 
further reduced with the IRAF and STSDAS routines. The images were all corrected for 
the intensity- and position-dependent charge transfer efficiency (CTE) by applying a linear 
ramp with a correction factor chosen according to the average counts of the sky background 
(Holtzman et al. 1995). The images taken with the same filter were then combined to remove 
cosmic rays and to produce a total-exposure map. 

To extract the Ha and [S il] line fluxes from the emission line images, it is necessary to 
subtract the continuum emission. As the giant Hn region NGC 5471 is very luminous, its 
Ha emission contributes significantly to the broad R-band (F675W) image. Therefore, we 
first subtracted a scaled Ha image from the F675W image, then scaled the resultant image 
and subtracted it from the line images. The scaling took into account the exposure times 
and the filter transmission curves (Biretta et al. 1996). 

The continuum-subtracted Ha and [S n] images were processed with the standard pro- 
cedures (Chen et al. 2000) to extract line fluxes and to produce the [S ll]/Ha ratio map. Due 
to the ~ 300 km s -1 redshift of M101, two corrections have been applied to the extracted Ha 
flux. First, as the filter transmission of the red-shifted Ha line is ~ 93% of the peak trans- 
mission, the extracted Ha flux is multiplied by a correction factor of 1.07. Second, the [Nil] 
A6548 line is red-shifted into the Ha bandpass at ~ 91% of the peak transmission. Adopting 
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the de-reddened [Nil] A6548+6584/H/3 ratio of 0.23 and the logarithmic extinction at H/3, 
c(H/3), of 0.23 for NGC5471B reported by Kennicutt & Garnett (1996), we find that the 
[N il] A6548 line contributes to ~ 2% of the total transmitted flux; therefore, an additional 
correction factor of 0.98 is applied to the extracted Ha flux. 

To suppress artificial fluctuations in the [S ll]/Ha ratio in regions where the nebular 
emission is faint (and hence noisy), we have clipped pixels with fluxes less than 3a above the 
sky background for both the Ha and [S il] images of NGC5471. The Ha, F675W, F547M, 
continuum-subtracted Ha, continuum-subtracted [S il], and [S ll]/Ha images of NGC5471 
are presented in Figure 1. The A-E components of NGC 5471 defined by Skillman (1985) are 
marked in Fig. la, and the [S ll]-bright shells are marked in Fig. If. Corresponding close-up 
images of NGC 5471B are shown in Figure 2. 

Stellar photometry was carried out using the APPHOT package in IRAF for the F547M 
and F675W images of NGC5471B. Because of the complexity in the stellar field, we manu- 
ally selected stellar sources in NGC 5471B and its surroundings and derived their apparent 
magnitudes, m F547M and m F675W . The uncertainties in the photometric measurements for- 
mally given by APPHOT are 0.1-0.2 mag. In most cases the true errors must be larger, if 
the uncertainties in aperture correction and sky background are included. The photometric 
uncertainties are generally larger in rriF675W than in m,F547M, because the nebular emission 
produces rapid variations in the background in the F675W image. Varying the background 
aperture may result in changes in mp675W by ~ 0.2 mag for most stars, and up to ~ 0.5 mag 
for stars near bright nebulosities. To reduce the errors introduced by the nebular contamina- 
tion in the F675W band, we have also carried out stellar photometry for the Ha-subtracted 
F675 W image, and designated this apparent magnitude as m,F675W'- Table 2 lists the re- 
sulting magnitudes and colors of the stellar sources in NGC 5471B. The stellar sources are 
marked on the F547M image in Figure 3a, and the Ha-subtracted F657W image is shown 
in Fig. 3b. 



2.2. KPNO 4 m Echelle Observations 

High-dispersion spectra of NGC 5471 were obtained with the echelle spectrograph on 
the KPNO 4 m telescope on 1999 June 30 and July I. The observations were made with the 
79-63 echelle grating and the 226-1 cross disperser. The 2048 x 2048 T2KB CCD detector 
was used to record the spectra. The pixel size is 24 /zm. To reduce the readout noise, the 
data were read out with the pixels binned by a factor of 2 in the spatial direction. The final 
data array thus has a scale of C//57 pixel -1 along the slit, and ~ 0.08A pixel -1 along the 
dispersion. 
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Two slit positions were observed, one centered on the hypernova remnant candidate 
NGC5471B (= Shell 1 in Fig. If) and the other on a [S ll]-bright shell to its north (Shell 2 
in Fig. If). For both observations, the slit was E-W oriented with a slitwidth large enough 
to cover the entire shell, 2"0 for NGC5471B and 1"5 for Shell 2. NGC5471B was observed 
with a total exposure time of 30 min, while Shell 2 only 10 min. The exact slit positions, 
shown in Fig. 2a, were recovered by matching the surface brightness profiles along the slits 
to those extracted from the WFPC2 Ha image. The accuracy of the position matching is 
better than / .'5. The instrumental FWHM, determined from Gaussian fits to the unresolved 
telluric lines, is 16.8 ± 0.5 km s" 1 for NGC5471B, and 12.8 ± 0.5 km s" 1 for Shell 2. 

The Ha+[Nn] echellograms of NGC5471B and Shell 2 and the [S il] AA6716, 6731 
echellogram of NGC 5471B are displayed in Figure 4. These new echelle data reveal high- 
velocity gas over a larger velocity range than that reported by Chu & Kennicutt (1986) 
because of an improvement in the S/N of data. 

The IRAF software was used for data reduction and analysis. We followed the same 
reduction procedures as described in Paper I. The sky-subtraction is more difficult for 
NGC5471B, because it is embedded in the giant Hn region NGC 5471, and the nebular 
emission fills the entire length of the slit. We have thus adopted the sky spectrum ex- 
tracted from the observations of MF83 reported in Paper I to remove the telluric lines 
in the spectrum of NGC5471B. The results are satisfactory, because MF83 is sufficiently 
close to NGC 5471B and the observations of MF 83 were made immediately after those of 
NGC 5471B. The sky-subtracted Ha line profile of NGC 5471B extracted over 2"9 along the 
slit is Hanning smoothed over five adjacent pixels and displayed in Figure 5. 



3. Supernova Remnant Candidates in NGC 5471 

SNRs are commonly diagnosed by their high [S ll]/Ha ratios, as S + can be collisionally 
excited in the cooling region behind SNR shocks and emit strongly in the [S n] AA6716, 
6731 lines (Raymond 1979). Although the [S ll]/Ha ratios of known SNRs are frequently 
observed to be 0.5-1.0 (Fesen, Blair, & Kirshner 1985), optical SNR surveys have used a 
lower threshold, 0.4-0.45, to identify SNR candidates because nebular backgrounds tend to 
lower the apparent [S ll]/Ha ratios (e.g., Long et al. 1990; Matonick & Fesen 1997). 

The [S ll]/Hci! ratio map of NGC 5471 (Fig. If) shows three shells with enhanced [S n] 
emission. The positions, sizes, and [S irj/Ha: ratios of these three shells are listed in Table 3. 
The apparent [S Iij/Ha ratios of these shells, 0.2-0.4, are higher than those of the background 
Hn region, ~ 0.1, but lower than the SNR identification threshold. These low observed 
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ratios are not surprising, given that the background Hn region is extremely bright. We have 
measured the average background emission from regions around the shells and subtracted 
the background from the [S n] and Ha images. The background-subtracted [S ll]/Ha ratios 
of the [S ll]-bright shells, also given in Table 3, are in the range 0.5-0.8, fully consistent with 
those of SNRs. 

The [S ll]-bright shells in NGC5471 have sizes of 50-80 pc, for a distance of 7.2 Mpc 
to M101 (Stetson et al. 1998). Interstellar shells with these large sizes and high [S ll]/Ha 
ratios can be either classical SNRs or superbubbles (Lasker 1977). To distinguish between 
these two possibilities, the stellar content and nebular dynamics need to be examined. 3 A 
shell is likely to be a superbubble, if it encompasses a high concentration of blue stars and 
the shell has an expansion velocity <C 100 km s _1 . Conversely, the existence of SNRs can 
be inferred, if no concentrations of blue stars are present or the shell expansion velocity is 
high, > 100 km s -1 . 

Shell 1, the brightest of the three [S ll]-enhanced shells in NGC5471, corresponds to 
the SNR NGC5471B. The HST WFPC2 images in Fig. 2 show concentrations of blue stars 
within the boundary of Shell 1, making the existence of a superbubble a distinct possibility. 
On the other hand, the echelle spectrum of Shell 1 (Fig. 5) shows high-velocity gas over a 
velocity range of ~ 620 km s _1 , which has been seen only in SNRs. It is thus most likely that 
Shell 1 contains a SNR but the physical conditions are complicated by the complex stellar 
environment. A detailed analysis of Shell 1 will be presented in §4. 

Shell 2 is at ~ 2" north of NGC5471B. With a size of 52 x 50 pc, it is the smallest of 
the three shells. It exhibits a clear ring morphology in both the Ha and [S n] lines. The 
F54 7M continuum image (Fig. 2c) shows sparse distribution of stars in the vicinity of Shell 2 
with no specific concentrations within the shell interior. The echellogram of Shell 2 (Fig. 4c) 
shows a very broad Ha line with high-velocity gas detected over a velocity range of ~ 450 
km s _1 . The lack of prominent concentrations of stars and the large expansion velocity in 
the high-velocity gas both support the existence of a SNR. Even if Shell 2 is a superbubble 
formed by a small OB association, it must have been recently accelerated by a SNR. 

Shell 3 is located at the low-density outskirts of NGC5471. Its size, ~ 70 pc, is com- 
parable to that of Shell 1. It has the highest [S ll]/Ha ratios among the three [S II]-bright 
shells, and its background-subtracted [S ll]/Ha ratio is within the higher range of SNRs. 
The F547M continuum image (Fig. lc) does not detect any bright stars or concentrations of 



3 The conventional SNR diagnostics, bright X-ray emission and nonthermal radio emission, are ineffective 
for SNRs in bright giant H n regions such as NGC 5471 because shocked fast stellar winds can produce X-ray 
emission and the thermal emission from the H n region overwhelms the nonthermal emission from SNRs. 
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stars in the interior or along the rim of Shell 3. It is unlikely that Shell 3 is a superbubble. 
Although no kinematic information is available for Shell 3 and the size of Shell 3 is larger 
than the average size of Galactic SNRs, we consider that Shell 3 is probably a bona fide SNR 
and its large size is a direct consequence of the low density in its interstellar environment. 

4. Analysis of Physical Properties of NGC 5471B 

4.1. Interstellar and Stellar Environment of NGC 5471B 

As shown in Fig. 1, NGC5471B has a complex environment. It is not trivial to separate 
the SNR from the background. Different morphologies are seen in the Ha and [S n] images. 
Since SNR shocks can be diagnosed by high [S ll]/Ha ratios, the [S n] image should illus- 
trate more clearly the boundary of the SNR. Indeed, the continuum-subtracted [S n] image 
provides the most well-defined limb-brightened shell morphology of NGC5471B. The shell 
is elongated with a size of 2'.'2 x 1"7, or 77 x 60 pc. A bright emission patch is present in the 
southwest corner. Along this direction, [S ll]-bright emission is still present exterior to the 
shell, making the SNR boundary somewhat ambiguous. Owing to this patch of [S ll]-bright 
emission, the [S ll]/Ha ratio map in Fig. 2f shows high ratios over a region ~ 50% larger 
than the [S n] shell. We will assume that the limb-brightened [S il] shell represents the SNR 
in NGC5471B. 

The continuum-subtracted Ha image shows an oval disk of emission with a slight limb- 
brightening at the west rim and a bright compact emission region at the southeast corner. 
This bright compact region, showing bright stellar emission in the continuum bands and 
having the lowest [S ll]/Ha ratio, is most likely a high-excitation Hll region. The oval disk 
of Ha emission, roughly coincident with the [S n] shell, is probably associated with the SNR. 
The SNR is superposed on a bright H II background; the average surface brightness at the 
main body of the SNR is only 2-4 times as bright as those of the surrounding regions. 

The F547M image shows many concentrations of stars in and around NGC5471B. 
Within the main body of NGC 5471B, two peaks of stellar emission are present at the eastern 
end. The brighter stellar source is coincident with the aforementioned bright high-excitation 
H II region; the fainter stellar source is not associated with prominent Ha emission, but the 
[S ll]/Ha ratio in its vicinity is low, indicating that it is also inside an Hll region. There is 
diffuse continuum emission along and interior to the south rim of the main body. Since the 
F5^7Mband does not contain any bright nebular emission lines and the distribution of this 
continuum light is different from that of the Ha emission, this diffuse continuum emission 
most likely originates from an extended distribution of unresolved stars similar to the "star 
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clouds" defined by Lucke & Hodge (1970) for OB associations in the Large Magellanic Cloud 
(LMC). 

Many peaks of stellar continuum emission are present in the vicinity of NGC5471B 
(see Fig. 3). Their brightnesses in the F547M, F675W, and Ha-subtracted F675W bands 
(m>F547M, rnF675W, and rriF675W') are given in Table 2. To determine the nature of these stellar 
sources, we present in Figure 6 a color-magnitude diagram of these sources using tufs^m 
and m,F675W'- To facilitate comparisons with normal stars, we have synthesized mps^M 
and rriF675W f° r stars of different spectral types and luminosity classes using the Kurucz 
(1993) stellar atmosphere models and the WFPC2 filter transmission curves (Biretta et al. 
1996). The synthetic photometry was calculated for a metallicity of 1/10 Z Q , a reddening 
of E(B — V) — 0.16 mag (Kennicutt & Garnett 1996), and a distance modulus of 29.3 
mag (= 7.2 Mpc, Stetson et al. 1998). In Fig. 6, the synthetic main sequence from 03 V 
downward is plotted in a thick line and the supergiants in a thin line with spectral types 
marked. It is immediately clear from the comparison that all identified stellar sources are 
either supergiants or composite with multiple stars. 

While the brightest stellar source, #18, is outside and to the south of the shell in 
NGC5471B, the next two brightest sources, #21 and #22, are the two stellar concentrations 
inside the shell. All these three bright sources have blue colors, and must be OB associations 
or clusters. Source #18, not surrounded by a recognizable Hn region, must be the oldest 
among these three sources. It is thus apparent that the SNR shell in NGC 5471B encompasses 
two young OB associations. 

4.2. Luminosity and Dynamics of the SNR Shell in NGC 5471B 

Several physical properties of the SNR shell in NGC5471B can be derived from the 
multi-order echelle observations. First of all, the expansion of the shell can be determined 
from the Ha line profile. The echelle image of the Ha and [Nil] lines of NGC5471B (Fig. 4) 
shows narrow cores and extended wings superposed on a faint stellar continuum. The pres- 
ence of the extended wings in the [N il] lines indicates that this component originates from 
shocked gas, as opposed to stellar atmospheric emission from Wolf-Rayet (WR) stars or 
luminous blue variables (LBVs). The Ha velocity profile of NGC5471B (Fig. 5), extracted 
over 2 / .'9 along the 2"-wide slit, can be fitted by a narrow component of FWHM = 41 ± 2 
km s" 1 and a broad component of FWHM = 148 ± 5 km s -1 . The widths of the narrow and 
broad components are in the ranges typical for Hn regions and SNRs, respectively (Chu & 
Kennicutt 1988, 1994b). Therefore, we assign the narrow component to the background Hn 
region and the broad component to the SNR. 
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The decomposition of SNR and Hn region components within the echelle aperture 
allows us to determine the Ha luminosity of the SNR from the HST WFPC2 Ha image of 
NGC5471B. The total Ha flux within the 2" x 2"9 aperture centered on the SNR shell is 
(2.6 ± 0.2) x 10 -13 ergs cm~ 2 s _1 . The error originates mostly from the imprecise aperture 
of the echelle observation which is affected by the seeing. The best spectral fit to the Ha 
velocity profile indicates that the broad component is 1.4 times as bright as the narrow 
component. Therefore, the Ha flux from the SNR is (1.5 ±0.1) x 10~ 13 ergs cm~ 2 s _1 . This 
flux is consistent with the value 1.5 x 10~ 13 ergs cm -2 s -1 estimated by Chu & Kennicutt 
(1986). For a distance of 7.2 Mpc and a c(H/?) of 0.23 (Kennicutt & Garnett 1996), the Ha 
luminosity of the SNR is (1.4 ± 0.1) x 10 39 ergs s _1 . 

The mass of the SNR can be derived from its Ha luminosity and density. Our echelle 
observation detected the [S n] AA6716, 6731 doublet. The velocity profiles of the [S il] 
doublets are likewise fitted by narrow and broad components. The [S n] A6716/A6731 ratio 
is 1.4 and 1.0 for the narrow and broad components, respectively. The [S n] doublet ratio 
for the narrow component is at the low-density limit, implying that the electron density of 
the Hn region is <C 50 cm~ 3 . The [S n] doublet ratio for the broad component, on the other 
hand, corresponds to an electron density of ~ 700 cm -3 , adopting an electron temperature 
of 13,000 K (Skillman 1985). This density and the above-derived Ha luminosity infer a mass 
of 4600 ±500 M in the SNR shell. This mass is consistent with but more accurate than the 
value 6500 ±3000 M estimated by Chu & Kennicutt (1986) because of the updated atomic 
constants and more direct measurement of the electron density used in our current work. 

Finally, the kinetic energy of the SNR shell can be estimated from the mass and ex- 
pansion velocity. The expansion velocity of the SNR cannot be measured directly because 
the SNR is not spatially resolved in the echelle data. As the integrated velocity profile of 
the SNR is weighted by the surface brightness of each emission region within the SNR, the 
FWHM of the integrated profile cannot be used to determine the expansion velocity of the 
SNR. Instead, we adopt the observed extreme velocity offsets of the broad line as a lower 
limit of the expansion velocity of the SNR. The approaching and receding sides of the SNR 
are detected up to about —290 km s -1 and ±330 km s _1 with respect to the systemic veloc- 
ity of the Hn region, respectively (see Fig. 5). Thus we adopt an expansion velocity of 330 
km s _1 for the SNR, and the kinetic energy of the SNR is 5.0 x 10 51 ergs. The kinetic energy 
estimated by Chu & Kennicutt (1986) is 20 times smaller because they used too small an 
expansion velocity for the SNR. 
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4.3. X-ray Emission and Hot Ionized Gas of NGC 5471B 

X-ray emission from NGC 5471 has been detected by ROSAT observations with both 
the Position Sensitive Proportional Counter (PSPC) and the High- Resolution Imager (HRI). 
The angular resolutions of the PSPC and HRI were not adequate to distinguish between 
point sources and diffuse emission from NGC 5471; however, the spectral shape of the PSPC 
spectrum of NGC 5471 is consistent with that of thin plasma emission (Williams & Chu 
1995). Using an ultradeep, 227 ks, ROSAT HRI image of M101, Wang, Immler, & Pietsch 
(1999) detected NGC 5471 at a count rate of 1.22 counts ks -1 , corresponding to an X-ray 
luminosity of 3 x 10 38 ergs s -1 (Wang 1999). 

We have plotted the X-ray contours extracted from the 227 ks ROSAT HRI image over 
the HST WFPC2 [S n] image of NGC 5471 in Figure 7. The X-ray contours are centered near 
NGC5471B. Given the combined pointing errors of the HST and ROSAT, a few arcsec, it 
is very likely that the X-ray source is centered at NGC5471B. The high-resolution Chandra 
observation of NGC 5471 (scheduled for Cycle 3) will be able to determine unambiguously 
the position and spatial extent of this X-ray emission. 

Wang (1999) used the observed X-ray emission from NGC 5471 in conjunction with the 
Sedov solution for a SNR, and computed a thermal energy of ~ 10 52 ergs for the hot gas in 
the NGC 5471B SNR. This amount of thermal energy is about twice as high as the kinetic 
energy we have derived for the 10 4 K ionized expanding shell of this SNR. 

5. Discussion 

The term "hypernova" has been vaguely used in theoretical contexts until recently 
when the connection between GRBs and supernovae was illustrated observationally through 
SN1988bw (Kulkarni et al. 1998) and SN 1997cy (Germany et al. 2000). Analyses of the 
light curves and spectral variations of these two supernovae indicate that the kinetic energy 
released in the supernova explosion was 3 — 5 x 10 52 ergs (Nakamura et al. 2001; Turatto et al. 
2000). Such explosion energies are more than an order of magnitude higher than the canon- 
ical supernova explosion energy of 10 51 ergs, qualifying these supernovae as "hypernovae" . 
Hypernovae without GRB counterparts have also been reported, e.g., SN1997ef (Mazzali, 
Iwamoto, & Nomoto 2000). Thus, we consider a supernova with an explosion energy greater 
than 10 52 ergs as a hypernova, and a SNR that requires such a high explosion energy as a 
hypernova remnant. 

The SNR in NGC5471B has been suggested to be a hypernova remnant because its 
X-ray emission implies a thermal energy that requires a supernova explosion energy greater 
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than 10 52 ergs (Wang 1999). Our high-resolution imaging and spectroscopic observations at 
optical wavelengths have allowed us to (1) resolve the shell structure of the SNR, (2) examine 
the stellar and interstellar environments of the SNR, (3) study the expansion dynamics of 
the SNR, (4) separate the Ha flux of the SNR from the background Hll region emission, 
and (5) determine the mass and kinetic energy of the SNR shell. Using these results we are 
now able to critically examine the nature of the hypernova remnant candidate NGC5471B. 

5.1. Explosion Energy Requirement of the SNR Shell in NGC 5471B 

NGC5471B possesses typical signatures of SNRs, e.g., nonthermal radio emission, high 
[S Il]/Hoj ratio, and X-ray spectrum consistent with those of thin plasma emission (Skillman 
1985; Williams & Chu 1995). It is most certain that NGC5471B contains a SNR. The HST 
WFPC2 images of NGC 5471B show a dominant, [S ll]-bright 77 x 60 pc shell, which is 
responsible for the enhanced [S Iij/Ha ratio. The echelle observations of NGC5471B have 
established that there is a broad velocity component associated with high-velocity shocks. 
From the comparison between the [S n] and Ha velocity profiles of NGC5471B (Fig. 4), it 
is clear that the broad velocity component is responsible for the [S ll]-enhancement. There- 
fore, the [S ll]-bright shell shown in the WFPC2 images must produce the broad velocity 
component in the echelle spectra. 

We can further confirm that the [S Il]-bright shell corresponds to the broad velocity 
component using a consideration of the brightnesses of the shell and the background Hll 
region. The Ha velocity profile shows that the broad component contributes to 1.4/(1+1.4) 
= 58% of the total Ha flux within the 2" x 2'.'9 echelle aperture. Using the WFPC2 Ha 
image and assuming that the shell is superposed on an Hll region with an average surface 
brightness similar to that of the surrounding regions, we find that the shell contributes to 
~ 54% of the total Ha flux within a 2" x 2'.'9 aperture centered on the shell. The similarity 
in the shell's and the broad component's contributions to the total flux within the same 
aperture strengthens the conclusion that the [S ll]-bright shell gives rise to the broad velocity 
component. The shell expansion velocity, ~ 330 km s~ x , is higher than those of most SNRs in 
the LMC (Chu & Kennicutt 1988), and certainly can produce the radio and X-ray signatures 
of SNRs by shocks 4 . Therefore, we believe that all observed SNR signatures originate from 
this [S ll]-bright shell, and will identify this shell as the SNR. 



4 Note that the shock velocities are higher than the observed expansion velocity of the optical SNR shell. 
For example, the velocity of an adiabatic shock is 4/3 times the shell expansion velocity. Furthermore, the 
optically detected material is likely associated with SNR shocks propagating through dense clouds, while the 
shocks in the intercloud medium are faster but their post-shock material is too tenuous to be detected. 
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As we have derived in §4.2, the expanding SNR shell in NGC5471B has a mass of 
4,600 M and a kinetic energy of 5.0 x 10 51 ergs. The kinetic energy in a mature SNR 
shell is much lower than 30% of the supernova explosion energy (Chevalier 1974); thus the 
observed kinetic energy in the SNR shell in NGC 5471B requires an explosion energy greater 
than 10 52 ergs. Moreover, Wang (1999) assumed that the X-ray emission from NGC5471B 
originates from the hot shocked gas in a SNR, and derived a thermal energy of ~ 10 52 ergs 
for the SNR. Both the optical and X-ray considerations require a supernova explosion energy 
> 10 52 ergs; therefore, the [S ll]-bright expanding shell in NGC5471B must be a hypernova 
remnant or energized by a burst of multiple supernova explosions. 



5.2. Comparison with Other SNRs in Nearby Star Formation Regions 

Was the [S ll]-bright shell in NGC 5471B produced by a single hypernova explosion? It 
is almost impossible to answer this question, especially since the shell is in an active star 
formation region with two prominent OB associations (see Figs. 3 & 6) projected within 
the shell, as described in §4.1. It is nevertheless possible to compare NGC5471B to other 
SNRs in nearby active star formation regions in order to determine whether the shell in 
NGC 5471 B is a common entity. We will use the sample of star formation regions in the 
LMC and NGC 604 in M33 for comparison. 

Convincing SNR signatures have been detected in at least 13 OB associations in the 
LMC (Chu 1997). Of these, eight contain SNRs showing signatures at all wavelengths, and 
all of them are in the outskirts of their associated OB associations. The other five contain 
SNRs within their superbubbles and show only X-ray signatures. These superbubbles have 
sizes 50-100 pc, comparable to that of the shell in NGC5471B; however, their expansion 
velocities are <C 100 km s -1 , much slower than that of the NGC5471B shell. The shell in 
NGC5471B is much more energetic than the normal superbubbles around OB associations. 

We may further dismiss the two prominent OB associations projected within the shell 
of NGC5471B as the providers of multiple supernovae. The existence of bright compact 
H II regions around these OB associations indicate that they are young and the surrounding 
interstellar gas has not been dispersed by SNR shocks. It is thus unlikely that the two 
OB associations in the NGC5471B shell energized the SNR shell. Moreover, the size and 
expansion velocity of the SNR shell implies an upper limit of 1 x 10 5 yr for the dynamic 
age. In order to have more than 10 normal supernovae from an OB association within such 
a small time span, these supernovae must have the shortest-lived progenitors, i.e., the most 
massive stars. Assuming that all 10 progenitors were 03 V stars with initial masses of 100- 
120 M Q and adopting Salpeter's (1955) slope for the initial mass function, we expect the 
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OB associations to contain 250 stars in the mass range of 25-100 M Q , corresponding to 
the later-type O stars (Schaerer & de Koter 1997). The mps^M magnitudes of the two OB 
associations projected within the NGC 5471B shell are equivalent to 70 and 45 06 V stars for 
sources #21 and #22, respectively. The total number of O stars in these two OB associations 
is most likely smaller than 115 because some of the stars may be supergiants that are several 
times brighter than their main-sequence counterparts. Therefore, the expected number of 
O stars exceeds the total number of O stars implied by m F5j , 7M by more than a factor of 2. 
Thus, the observed OB associations in the shell of NGC 5471B are unable to produce enough 
supernovae within the dynamic age and power the shell. 

It would be appropriate to compare NGC5471B to the giant Hn region 30 Doradus in 
the LMC, as 30 Dor is the most active star formation and has the highest concentration of 
massive stars in the LMC. Figure 8 displays NGC5471B and 30 Dor at the same physical 
scale over an area of 200 pc x 200 pc, adopting a distance of 50 kpc to the LMC (Feast 
1999). Within this region, the Ha luminosity of NGC 5471B is 4.6 x 10 39 ergs s _1 , and 30 Dor 
6.2 x 10 39 ergs s -1 (Kennicutt & Hodge 1986; Mathis, Chu, & Peterson 1985). The similarity 
in the Ha luminosity implies a similarity in the ionizing power of stars encompassed in these 
two regions. It would then be truly a fair comparison between NGC 5471B and 30 Dor. 

Despite the similarity in ionizing power, the nebular structure of NGC5471B is very 
different from 30 Dor. NGC 5471B is dominated by the SNR shell, while 30 Dor is dominated 
by several shell structures of comparable sizes. The kinematic structure of 30 Dor has been 
mapped in great detail by Chu & Kennicutt (1994a). The central shell around the R136 
cluster in 30 Dor is much smaller (30 pc in diameter) and expands much more slowly (~ 130 
km s _1 ) than the shell in NGC5471B. The five largest shells in 30 Dor, designated as Shells 
1-5 by Wang & Helfand (1991), show a common kinematic structure: a slowly expanding 
shell with V exp <C 100 km s _1 superposed by additional shocked clouds with velocity offsets 
up to ±300 km s _1 . The Ha fluxes from these shells in 30 Dor are dominated by the slow 
expanding shell component. We use Shell 1 in 30 Dor as an example. The echelle observations 
along four slit positions in 30 Dor-Shell 1, marked in Fig. 8, are used to simulate an integrated 
Ha velocity profile of 30 Dor-Shell 1, shown in Figure 9. The "pseudo-integrated" profile 
of 30 Dor-Shell 1 shows a narrow core and broad wings, but the broad wings are highly 
asymmetric and the broad component contributes to only 22% of the total flux. 

The fast expanding shell in NGC 5471B is unmatched by any comparable shells in the 
LMC. The giant Hn region NGC 604 in M33 has also been studied in detail with similar 
echelle observations (Yang et al. 1996). Its kinematic structure is very similar to that of 
30 Dor, and it does not have any shell with an expansion velocity or luminosity matching 
those of the shell in NGC5471B. Therefore, it becomes clear that the fast expanding shell in 
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NGC5471B is truly outstanding and unique. It could not have been produced by multiple 
supernovae as in the giant Hn regions 30 Dor and NGC604. Therefore, we conclude that it 
is very likely that the energetic SNR in NGC 5471B was produced by a hypernova. 

This research is supported by the NASA grant STI6829.01-95A. 
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Fig. 1.— HST WFPC2 images of NGC5471 in (a) Ha, (b) F675W, (c) F547M, (d) 
continuum-subtracted Ha, (e) continuum-subtracted [S n], and (f) [S ll]/Ha. The grey-scale 
images in (a)-(e) are displayed with square- root image transfer functions. The [S ll]/Ha 
ratio map in (f) is displayed with a linear image transfer function, as shown in the greyscale 
bar to its right. 
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Fig. 2.— HST WFPC2 images of NGC5471B in (a) Ha, (b) F675W, (c) F547M, (d) 
continuum-subtracted Ha, (e) continuum-subtracted [S il], and (f) [S ll]/Ha. The grey-scale 
images in (a)-(e) are displayed with square- root image transfer functions. The [S ll]/Ha 
ratio map in (f) is displayed with a linear image transfer function, as shown in the greyscale 
bar to its right. The lines in (a) mark the widths and locations of the E-W oriented slit of 
the echelle observations. 
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Fig. 3.— HST WFPC2 images of NGC5471B in (a) F547M and (b) Ha-subtracted F657W. 
Stellar sources are identified and marked in (a) 
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Fig. 4.— KPNO 4 m echellograms of (a) the He*+[Nn] A6583 lines, and (b) the [S n] AA6716, 
6731 lines along an E-W slit centered on NGC5471B (=Shell 1), and (c) the Ha+[Nn] A6583 
lines along an E-W slit centered on Shell 2. The C-component of NGC5471, to the east of 
NGC5471B, is also detected and marked in (a) and (b). The bright source to the east of 
Shell 2 is a bright compact H II region to the north of NGC 5471C. The images are displayed 
with square-root image transfer functions, as indicated by the greyscale bars on the right. 
The telluric lines marked by © are, from left to right, Ha A6562.8, OH A6568.8, and OH 
AA6577.2, 6577.4 lines. Both NGC5471B and Shell 2 have an average heliocentric velocity 
of 297 ±3 km s -1 . 
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Fig. 5. — Ha velocity profile of NGC 5471B extracted from a 2" x2'.'9 echelle aperture centered 
on NGC5471B. The profile (in solid line) is fitted by two Gaussian components (in dotted 
and dashed lines). The residual of the fit, displayed in the lower panel, shows that high- 
velocity gas is detected over more than 600 km s _1 . The relative velocity in the horizontal 
axis is referenced to the heliocentric velocity of NGC5471B, 297 km s _1 . 
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Fig. 6. — Color-magnitude diagram of the stellar sources in the vicinity of NGC5471B. The 
stellar sources in Table 2 are plotted as •. The Ha-corrected F675W magnitude, rriF675W', 
is used to reduce photometric errors caused by Ha contamination in the F675W band. Loci 
of the synthetic main sequence from 03 V downward are plotted in a thick line and the 
supergiants in a thin line, for a metallicity of 1/10 Z@, a reddening of E(B — V) — 0.16 
mag, and a distance modulus of 29.3 mag. The spectral types of the synthetic supergiants 
are marked. 
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Fig. 7.— (a) ROSAT HRI image of NGC5471. This HRI image is extracted from a 227 ks 
observation, and has been smoothed with a Gaussian of a = 3". (b) HST WFPC2 [S n] 
image overlaid with the HRI X-ray contours at levels of 25%, 50%, 75% and 90% of the peak 
value. 




Fig. 8.— (a) HST WFPC2 Ha image of NGC 5471B. (b) Ground-based Ha image of 30 Dor 
binned to simulate a HST WFPC2 image of 30 dor at the distance of M101. Both images 
have a linear dimension of 200 pc x 200 pc, and are displayed with square-root image transfer 
functions. 30 Dor-Shell 1 and the echelle slit positions from Chu & Kennicutt (1994a) are 
marked over 30 Dor in (b). 
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Fig. 9. — Pseudo-integrated Ha velocity profile of Shell 1 in 30 Dor. The profile (in solid 
line) is fitted by three Gaussian components (in dotted, short-dashed, and long-dashed lines). 
The residual of the fit is displayed in the lower panel. The relative velocity in the horizontal 
axis is referenced to the heliocentric velocity of 30 Dor, 275 km s _1 . 
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Table 1: Journal of HST WFPC2 Observations of NGC5471 



Filter 


A a 
(A) 


AA a 

(A) 


Band 


Exposure 
(sec) 


Image Root Name 


F656N 


6562 


22 


Ha 


1 x 180 


u4dn0305r 










2 x 600 


u4dn0306r, u4dn0307r 


F673N 


6733 


47 


[S II] 


3 x 700 


u4dn0301r, u4dn0302r, u4dn0303r 


F547M 


5454 


487 


Stromgren y 


2 x 100 


u4dn030fr, u4dn030gr 










2 x 600 


u4dn030dr, u4dn030er 


F675W 


6997 


889 


WFPC2 R 


2 x 50 


u4dn030ar, u4dn030br 










2 x 400 


u4dn0308r, u4dn0309r 



a The mean wavelengths (A) and widths (AA) of the filters are taken from the WFPC2 Instrument Handbook 
(Biretta et al. 1996). 
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Table 2: Photometry of Stellar Sources in NGC5471B 



ID 


R.A. 


Decl. 


n^F547M 


™>F675W 


TTlF675W' a m F547M — mF675W> 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) (mag) 



1 


14 


04 


28.94 


54 


23 


53.9 


22.68 


± 


0.08 


24.11 


± 


0.25 


23.52 


± 


0.12 


-0.85 


± 


0.14 


2 


14 


04 


29.00 


54 


23 


53.0 


22.23 


± 


0.10 


23.15 


± 


0.20 


22.92 


± 


0.13 


-0.69 


± 


0.16 


3 


14 


04 


29.02 


54 


23 


53.6 


21.48 


± 


0.04 


22.36 


± 


0.05 


22.24 


± 


0.05 


-0.77 


± 


0.06 


4 


14 


04 


29.03 


54 


23 


50.9 


23.80 


± 


0.47 


23.58 


± 


0.56 


22.69 


± 


0.11 


1.11 


± 


0.48 


5 


14 


04 


29.04 


54 


23 


50.2 


22.99 


± 


0.21 


21.91 


± 


0.10 


22.64 


± 


0.10 


0.36 


± 


0.23 


6 


14 


04 


29.06 


54 


23 


50.7 


22.49 


± 


0.10 


22.49 


± 


0.15 


22.40 


± 


0.08 


0.09 


± 


0.13 


7 


14 


04 


29.07 


54 


23 


51.8 


20.86 


± 


0.03 


20.88 


± 


0.05 


21.07 


± 


0.03 


-0.21 


± 


0.05 


8 


14 


04 


29.09 


54 


23 


53.8 


21.97 


± 


0.06 


22.65 


± 


0.08 


22.78 


± 


0.10 


-0.81 


± 


0.12 


9 


14 


04 


29.10 


54 


23 


53.4 


21.80 


± 


0.07 


22.51 


± 


0.09 


22.39 


± 


0.08 


-0.59 


± 


0.10 


10 


14 


04 


29.10 


54 


23 


54.3 


23.08 


± 


0.19 


23.55 


± 


0.23 


24.01 


± 


0.28 


-0.93 


± 


0.34 


11 


14 


04 


29.10 


54 


23 


54.8 


22.90 


± 


0.09 


24.62 


± 


0.56 


23.82 


± 


0.17 


-0.92 


± 


0.20 


12 


14 


04 


29.12 


54 


23 


50.1 


23.53 


± 


0.31 


23.29 


± 


0.32 


23.99 


± 


0.36 


-0.46 


± 


0.48 


13 


14 


04 


29.12 


54 


23 


51.4 


22.75 


± 


0.18 


23.71 


± 


0.86 


23.23 


± 


0.36 


-0.48 


± 


0.41 


14 


14 


04 


29.15 


54 


23 


53.3 


21.47 


± 


0.05 


22.21 


± 


0.09 


22.24 


± 


0.08 


-0.77 


± 


0.09 


15 


14 


04 


29.16 


54 


23 


54.2 


22.00 


± 


0.05 


22.71 


± 


0.09 


22.86 


± 


0.08 


-0.86 


± 


0.09 


16 


14 


04 


29.18 


54 


23 


51.4 


23.97 


± 


0.81 


23.31 


± 


0.85 


23.97 


± 


0.83 


0.00 


± 


1.16 


17 


14 


04 


29.21 


54 


23 


53.1 


22.68 


± 


0.26 


20.94 


± 


0.08 


21.79 


± 


0.11 


0.89 


± 


0.28 


18 


14 


04 


29.27 


54 


23 


51.3 


19.73 


± 


0.01 


20.09 


± 


0.02 


20.16 


± 


0.01 


-0.43 


± 


0.01 


19 


14 


04 


29.28 


54 


23 


50.0 


23.22 


± 


0.11 


24.24 


± 


0.28 


24.47 


± 


0.25 


-1.25 


± 


0.27 


20 


14 


04 


29.28 


54 


23 


54.5 


22.49 


± 


0.06 


22.96 


± 


0.16 


23.26 


± 


0.13 


-0.77 


± 


0.14 


21 


14 


04 


29.29 


54 


23 


52.4 


20.14 


± 


0.02 


19.66 


± 


0.03 


20.63 


± 


0.04 


-0.49 


± 


0.05 


22 


14 


04 


29.29 


54 


23 


52.9 


20.66 


± 


0.04 


20.32 


± 


0.07 


21.07 


± 


0.07 


-0.41 


± 


0.08 


23 


14 


04 


29.33 


54 


23 


50.2 


22.69 


± 


0.07 


23.82 


± 


0.27 


23.72 


± 


0.14 


-1.03 


± 


0.16 


24 


14 


04 


29.34 


54 


23 


53.8 


23.11 


± 


0.17 


23.40 


± 


0.21 


23.79 


± 


0.20 


-0.68 


± 


0.26 


25 


14 


04 


29.34 


54 


23 


54.3 


22.93 


± 


0.10 


24.05 


± 


0.36 


24.06 


± 


0.23 


-1.13 


± 


0.25 


26 


14 


04 


29.39 


54 


23 


51.8 


21.10 


± 


0.03 


20.89 


± 


0.03 


21.32 


± 


0.03 


-0.22 


± 


0.04 


27 


14 


04 


29.44 


54 


23 


53.1 


21.31 


± 


0.03 


21.84 


± 


0.06 


21.97 


± 


0.04 


-0.66 


± 


0.05 


28 


14 


04 


29.46 


54 


23 


51.7 


23.70 


± 


0.19 


23.59 


± 


0.24 


24.54 


± 


0.38 


-0.84 


± 


0.43 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arminutes, and arcseconds. 
a Apparent magnitudes measured from the Ha-subtracted F675W image of NGC5471B. 
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Table 3: [S ll]-Bright Shells in NGC5471 



Shell 


R.A. 


Decl. 


Angular Size 


Linear Size a 


Apparent 


Corrected 13 


ID 


(J2000) 


(J2000) 


(") 


(PC) 


[S n]/Ra 


[S II] /Ha 


1 


14 04 29.27 


54 23 52.5 


2.2 x 1.7 


77 x 60 


0.33 ± 0.03 


0.5 ± 0.1 


2 


14 04 29.20 


54 23 54.7 


1.5 x 1.4 


52 x 50 


0.23 ± 0.03 


0.5 ± 0.2 


3 


14 04 30.08 


54 23 43.0 


2.0 x 1.9 


70 x 67 


0.40 ± 0.05 


0.8 ± 0.2 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arminutes, and arcseconds. 
a For a distance of 7.2 Mpc to M101 (Stetson ct al. 1998). 

fc These corrected ratios are derived with an average background H n region emission subtracted. 



